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1.   Introduction  and  discussion  of  results 

In  1844  Scott  Russell  [1]'""  reported  his  experimental 
observations  on  the  solitary  wave,  a  wave  consisting  of  a 
single  elevation  which  propagates  without  change  of  form. 
(The  wave  motion  is  in  two  dimensions  only.)   Later 
Boussinesq  [2]  and  Rayleigh  [3]  independently  gave  deri- 
vations of  the  approximate  form  and  velocity-amplitude 
relationship  for  such  a  wave  in  shallow  water.   Then  by  a 
slight  modification  of  Rayleigh' s  method,  Korteweg  and 
DeVries  [4]  obtained  periodic  waves  of  permanent  type, 
which  include  the  solitary  wave  as  a  special  case  when 
the  wave  length  becomes  infinite.   Gwyther  [5]  obtained 
Rayleigh' s  results  by  a  slightly  different  method,  and 
McCowan  [6]  found  more  accurate  results  by  guided  guess- 
ing.  The  question  of  the  existence  of  the  solitary  wave 
In  water  of  finite  depth  has  been  considered  by  Weinstein 
[7] ,  but  the  existence  of  such  a  wave  has  not  yet  been 
proved. 

All  the  results  mentioned  above  apply  to  irrota- 
tional,  two-dimensional  motion  of  an  incompressible,  ir.- 
viscid  fluid  over  «  horizontal  bottom,   Rayleigh' s  treat- 
ment involves  an  iteration  scheme  of  a  peculiar  kind  and 
leads  to  a  differential  equation  for  the  wave  profile . 


Numbers  in  square  brackets  refer  to  the  bibliography  at 
the  end  of  the  article. 


The  theory  given  by  Bousslnesq  involves  a  ntunber  of  phys- 
ical assumptions  in  addition  to  those  of  the  basic  hydro- 
dynamical  theory;  it  also  leads  to  a  differential  equa- 
tion for  the  wave  profile.   Both  of  these  methods  assume 
that  the  depth  of  the  water  is  small  compared  to  some 
horizontal  dimension,  and  they  might  be  interpreted  as 
developments  of  the  whole  problem  in  powers  of  the  ratio 
of  the  depth  to  some  horizontal  dimension,  such  as  wave 
length.   However,  because  these  procedures  are  so  unsys- 
tematic, it  is  not  clear  that  they  are  equivalent  to  such 
developments,  nor  to  what  order  of  approximation  the  so- 
lutions obtained  are  valid.   The  method  of  proceeding  to 
higher  -approximations  is  also  obscure. 

The  object  of  the  present  investigation  is  to  dis- 
cuss waves  of  permanent  type"  in  shallow  water  by  a 
method  in  which  the  character  of  the  approximation  is 
quite  clear,  and  which  is  capable  of  being  carried  out  to 
Include  terms  of  any  desired  order.   The  method  consists 
in  expanding  the  solution  of  the  exact  hydrodynamic  prob- 
lem systematically  in  powers  of  a  dimensionless  parameter 
cf  =  ojh  where  h  is  the  depth  of  the  undisturbed  fluid  and 
CO  is  the  curvature  at  some  point  on  the  surface.   The  ex- 
pansions are  inserted  into  the  equations  of  motion  and 
the  boundary  conditions,  and  coefficients  of  like  powers 
of  <3   are  equated.   The  variables  are  chosen  in  such  a  way 
that  the  terms  of  zero  order  in  the  expansion  in  powers 
of  a   satisfy  the  well-known  equations  of  the  nonlinear 
shallow  water  theory,  which  are  analogous  to  the  equa- 
tions of  gas  dynamics."" 


Actually  stationary  solutions  are  fo\ind;  progressive 
waves  may  be  obtained  from  them  by  adding  a  constant 
velocity  to  the  fluid. 

"These  equations  were  first  derived  in  this  way  by  K.  0. 
Priedrichs  [9],  to  whom  this  method  is  due. 


It  is  easily  shown  that  the  only  solutions  of  these 
equations  for  the  first  approximation  (satisfied  by  the 
terms  independent  of  a)    which  are  of  permanent  form  are 
the  constant  or  piecewise  constant  (shock  tyipe)    solu- 
tions.  However,  the  equations  for  the  second  approxima- 
tion (satisfied  by  the  coefficients  of  d)    have  solutions 
which  yield  periodic  waves  of  permanent  form,  similar  to 
those  of  Korteweg  and  DeVries,  and  also  solitary  waves, 
similar  to  those  of  Rayleigh  and  Boussiilesq. 

The  solution  of  permanent  form,  as  given  by  the 
first  and  second  approximations,  for  the  irrotational, 
two-dimensional  motion  of  an  incompressible,  inviscid 
fluid  of  mean  depth  h  over  a  horizontal  bottom  is 

(1.11    n   =  n    '      +(^    -  rj    '    )cn^[^  2Ft  (k)]  , 
^    '     *    *min    ^imax   imxn'         ^x        1^    '^' 

P  =  9Zin   -  7)  > 

V  =  0, 
_u_  ^  ^max  _   _  ,;?max  _  ^Wi,   2.x  ^^    ,,  ., 

■X         4  'Zmin  "^^^ 

In  these  equations  r^  is  the  surface  elevation,  measured 

up  from  the  bottom;  x  is  the  horizontal  coordinate  and  y 

is  the  vertical  distance  above  the  bottom;  n_  ^  and  7?_.^ 

'  max     '  min 

are  the  maximum  and  minimum  svirface  elevations,  respec- 
tively; A  is  the  wave  length,  p  the  pressure,  p  the  fluid 
density,  g  the  acceleration  of  gravity,  v  the  vertical 
velocity,  u  the  horizontal  velocity;  F-,  the  complete  el- 
liptic integral  of  the  first  kind  of  modulus  k;  en  the 
Jacobi  elliptic  fimction  of  modulus  k;  L  and  k  are  cer- 
tain parameters  given  by  the  last  two  equations  of  (5.24), 


This  solution  yields  a  two-parameter  family  of  peri- 
odic waves  as  does  that  of  Korteweg  and  DeVries .   The 
parameters  ^i  and  71    -      are  subject  only  to  the  condi- 

tions 

a.2)  o<2j^<i,     W>2.!i^, 

which  both  follow  from  the  definition  of  h  or  the  rela- 
tions among  n    ,  ^jj,j_j,j  L?  and  k.   For  any  values  of  y^ 

and  V^.^   satisfying  the  inequalities  above,  the  sxirface 
'mm       o      o  rL 

profile  is  periodic  in  x  with  the  wave  length  X.      The 

height  of  the  crest  above  the  mean  height  is  greater  than 

the  depth  of  the  trough  below  the  mesin  height.   The  crest 

is  also  narrower  than  the  trough.   The  wave  is  thus  not 

symmetric  about  the  mean  height  as  it  is  in  the  linear 

theory.   This  is  shown  in  Figure  1,  where  one  wave 

length  of  a  typical  wave  profile  is  plotted. 

Figure  2  shows  the  contour  lines  of  ^  as  a  function 

"^min     '?max 
of  — c—  and  — r—  .  From  this  figure  it  can  be  seen  that 

A  increases  as  either  j)    .      or  n     increases.   The  quali- 

'  mxn    *  max 

tative  behavior  of  this  function  is  the  same  as  that  of 
the  f -unction  obtained  by  Korteweg  and  DeVries,  but  quan- 
titatively it  is  slightly  different.   Thus,  although  for- 
mally the  equation  for  the  wave  profile  obtained  by 
Korteweg  and  DeVries  is  the  same  as  the  first  of  equa- 
tions (1.1),  the  two  px-oxiles  are  slighiily  differeiio  be- 
cause of  the  different  expressions  for  Mr? ^>V    •^)« 

tmax  'mm 

Anothci'  difference  between  their  solution  aiid  equationo 
(1.1)  occurs  in  the  expression  for  u,  which  they  find  to 
be  independent  of  both  x  and  y.   Prom  considerations  of 
conservation  of  mass,  that  is  much  less  reasonable  than 
the  dependence  on  x  which  is  given  by  the  foiirth  of  equa- 
tions (1.1).   The  second  of  equations  (1.1)  yields  the 
pressure  in  the  fluid,  which  is  seen  to  be  given  by  the 


hydrostatic  expression  even  in  the  second  approximation. 
This  expression  therefore  seems  to  be  very  accurate.   In 
the  Korteweg-DeVries  solution  the  pressure  is  not  ob- 
tained. 

Prom  Pig\ire  2  it  can  be  seen  that,  for  a  fixed  value 

of   ?•   ,  as   y^   decreases,  A  decreases  until  X  =  0  when 

is^  =  2  -  -2^  .  In  this  limiting  case  the  amplitude 
may  be  finite  but  the  wave  length  is  zero,  so  that  the 
surface  is  everywhere  discontinuous.   On  the  other  hand, 

if  -^5  increases,  indefinitely,  so  does  X.   While  this 

occurs,  the  troughs  become  v/ider  and  the  crests  higher 
and  relatively  narrow. 

It  was  hoped  that  the  solution  would  impose  an  upper 

limit  on  °^5^  ,  and  that  for  the  limiting  solution  the 

siirface  slope  would  be  discontinuous  at  the  crests. 
Then,  as  shown  by  Stokes  [8,  p.  418],  an  angle  of  120 
would  be  formed  at  the  crests.   However,  to  the  present 
order  of  approximation  these  results  are  not  obtained. 
Further  research  leading  to  these  results  would  certainly 
be  worth  while. 

y/hen  -2^  and  ^^   are  nearly  equal  to  one,  both  the 

Korteweg-DeVries  solution  and  equabions  (1.1)  reduce  to 
the  cosine  solution  of  the  linear  shallow  water  theory. 

Ysfhen  ^^   =  1  and  ^^   is  greater  than  one,  but  otherwise 

unrestricted,  the  wave  length  becomes  infinite  in  both 
solutions.   In  this  case  equations  (1.1)  become  (see  Fig- 
ure 3  for  wave  profile) 


P  =  PgC'?  -  7)' 
V  =  0, 


u 


v/gS 

\  =  CX)  . 

Equations  (1.3)  represent  a  solitary  v;ave .   The  profile, 
given  by  the  first  of  equations  (1.3),  is  the  same  as 
that  found  by  Boussinesq  [2] .   The  Korteweg-DeVries  solu- 
tion reduces,  in  this  case,  to  the  solitary  wave  found  by 
Raylelgh.   This  solution,  for  the  profile,  differs  from 
the  first  of  equations  (1.3)  only  by  having  the  additional 

factor  i/-r multiplying  the  argument  of  the  hyperbolic 

If  *  max 
secant.   For  waves  of  small  amplitude  this  factor  is 

nearly  one,  and  thus  in  this  case  the  Korteweg-DeVrles- 

Rayleigh  solution  for  the  profile  practically  agrees  with 

the  first  of  equations  (1.3).   Both  the  Korteweg-DeVries - 

Raylelgh  solution  and  the  Boussinesq  solution  yield  a 

horizontal  velocity  Independent  of  x   and  y.   As  mentioned 

above,  considerations  of  conservation  of  mass  indicate 

that  the  dependence  of  the  velocity  on  x  as  given  by  the 

fourth  of  equations  (1.3)  is  more  reasonable  than  the 

constant  velocity.   However,  for  x  infinite,  the  fourth 

^max  ,  T 

of   equations    (1.3)    yields  u  =  \/gh*  ,    and   the 

2 

velocity  Vgn   '  given  by  Korteweg-DeVries -Raylelgh  agrees 
with  this  to  first  order  in  the  relative  amplitude 
'Zmax  -  ^ 


Thus,  if  the  water  at  inTlnity  is  at  rest,  the  prop- 
agation speeds  of  the  wave,  as  given  by  the  two  solutions, 
agree  to  first  order  in  the  relative  amplitude  of  the  wave. 
Kov;ever,  according  to  the  fourth  of  equations  (1.3)  the 
water  under  the  crest  would  be  moving,  while  in  the 
Korteweg-DeVries-Rayleigh  solution  it  would  be  stationary. 

2 .   Derivation  of  the  shallow  water  theory 

In  this  section  the  exact  equations  of  hydrodynamics 
and  boundary  conditions  are  employed.   Viscosity  and  com- 
pressibility are  not  considered,  and  the  bottom,  given  by 
y  =  d(x,z),  is  assumed  to  be  rigid.   This  implies  that  the 
normal  component  of  velocity  is  zero  on  the  bottom.   The 
pressure  p(x,y,z)  is  assumed  to  be  a  constant  (which  is 
taicen  to  be  zero)  on  the  surface  y  =  ?^(x,z).   In  addition, 
the  motion  is  assumed  to  be  irrotational.   V/ith  these  as- 
siiRiptions,  and  the  kinematic  condition  that  a  particle  on 
the  surface  remains  on  the  surface,  the  equations  of  motion 
and  boiondary  conditions  are  as  follows  (the  y-axis  is  ver- 
tically upward;  u,  v,  w  are  the  x,  y,  z  velocity  compo- 
nents; p  is  the  density  and  g  the  acceleration  of  gravity): 


(2.1)  u+v+w„=0, 

^    '  X    y    z 


^t  "*"  ^^x  "^  ^^y  "^  ^z  ~  "^x^P' 

^t  "^  ^^x  ^  ^%  "^  ^z  ^  -Py/p-g' 

w,  +  uw_  +  vw,^  +  ww^  =  -p„/p, 

t     X  y  ^     ^ 


^y  =  ^z'    ^z  =  ^'^x'    ^x  =  ^y' 


^t  ^   ^^x  "^  ""^z   "  ^     at  y  =  T2, 
p  =  0     at  y  =  ?2» 


ud   -  V  +  wd„  =0     at  y  =  d. 
X         z 
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A  transformation  to  dimenslonless  variables  will  now 
be  performed.   The  constants  h  and  co,  with  dimensions  of 

length  and  reciprocal  length  respectively,  will  be  intro- 

2   2 
duced.   Later  the  dimenslonless  constant  cf  =  co  h  v/ill  be 

an  expansion  parameter,  and  it  Is  small  when  h,  the  typi- 
cal length  in  the  vertical  direction,  is  small  compared 
to  go"  ,  the  typical  length  in  the  horizontal  direction. 
For  this  reason  the  equations  obtained  when  cT  is  small 
are  called  the  equations  of  the  shallow  water  theory. 

The  transformation  equations,  defining  the  new  di- 
menslonless variables,  are 

(2.2)        X  =    ccco"-'-,         r  =   \/gEc:ot,        u  =  \/^  u, 

y  =  ph,  7  =  hY,  V  =  v'gh'(£Oh)"-'-v, 

z   =  Y^~    »  ^  =  hH>  w  =   \/gh!  w, 

2   2 
p   =  pghtr,  cf  =  CO  h    . 

It  is  essential  that  the  scale  factors  in  the  horizontal 
and  vertical  directions  are  different.   V/hen  these  ex- 
pressions are  used  in  equations  (2.1),  one  obtains 


(2.3)  du^  +  "^B  "^  ^^r  ^  °' 

<?[u^  +  uu^  +  wu^  +  rr^]  +  vu     =  0, 

d[v^    +   UV^    +   WV^  +   TT       +1]  +   vv      =0, 

cf[w^   +  uw„    +  WW.,  +  -rr.,]  +  vw„    =0,        "^^ 

cJ[Y^  +  uY^   +  wY^]    =  V        at  p    =  Y, 

■rr  =  0        at  p    =  Y, 

a[vE     +  wH  J    =  V        at  p    =  H. 
a  Y 


In  order  to  solve  these  equations,  it  will  be  as- 
sumed that  u,  V,  w,  TT,  and  Y  can  be  written  in  pov/er 
series  in  a.      These  series  will  then  be  inserted  into  the 
above  equations,  and  the  coefficients  of  like  pov/ers  of  cr 
equated.   The  coefficient  of  d   in  each  series  has  the 
superscript  n  (e.g.  u^)  and  is  a  function  of  a,  p,  y,  and 
Z,    except  the  coefficients  i     which  are  independent  of  p . 

Thus  inserting  the  series  and  retaining  terms  of 
zero  order  in  d   in  the  above  equations,  one  has 

(2.4)  Vp  =  0' 

.%;  =  0. 
.%;  =  0, 

o    o     o    o     o    o 
'^P  =  V'    \   =  ""cc'         ^0.  =  ^p' 

v°  =  0   at  p  =  Y°, 
^°  =  0   at  p  =  Y°, 


v°  =  0   at  p  =  H. 


These  equations  yield 


(2.5)  v°(a,p,Y,'C)  =  0, 


w°  =w°(fl6,Y,r), 


u°  =  \x°{oc,r,z) , 


Tr°(a,Y°,Y,r)  =  0, 
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Nov/,  retiirning  to  equations  (2.3)  and  equating  coeffi- 
cients of  the  first  pov/er  of  d  to  zero  one  obtains,  after 
using  equations  (2.5), 

(8.6)  <  +  "r  =  -"r 

o  o  o    o  o    o    ^ 

U^  +  U  U   +  W  U,  +  TT   =  0, 

t;      a,  Y    a. 

Tr°  .  1  =  0, 
o    o  o    o  o    o   „ 

W_  +  U  W   +  W  W  ,  +  TT.,  =  0 , 

Z  oc  X  Y 

Y°  +  u°Y°  +  w°Y°  =  v^   at  p  =  Y°, 
r     a.      Y  p     ' 

u°H^  +  w°H  =  v-'-   at  p  =  H. 
The  first  of  these  equations  can  be  integrated,  yielding 

(2.7)  v^  =  -(u°  +  w°)p  +  P(a,r) 

where  F{oc,z)    can  be  determined  from  the  last  of  equations 
(2.6).   Thus,  after  determining  P(a,r) , 

(2.8)  v^  =  -p(u°  +  w°)  +  (u°H)^  +  (w^H)^. 

Similarly  the  third  of  equations  (2.6)  can  be  inte- 
grated and  the  integration  constant  determined  from  the 
last  of  equations  (2.5).   This  gives  the  hydrostatic  ap- 
proximation 

(2.9)  iT°((X,p,Y,2;)  =  Y°(a,Y,2r)  -  p, 


which  is  talien  by  Lamb  [8,  p.  254]  as  the  starting  point 
of  the  shallow  water  theory.   Here  this  hydrostatic  pres- 
sure relation  is  automatically  satisfied  by  the  solution 


11 


in  the  first  approximation.   After  using  equations  (2.8) 
and  (2.9)  in  the  second,  foiirth,  and  fifth  of  equations 
(2.6) ,  one  has 

(2.10)  u°  +  u°u°  +  w°u°  +  Y°  =  0, 

r      a      Y    a 

w°  +  u°w°  +  w°w°  +  Y°  =  0, 

Y°  +  [u°(Y°  -  H)]^  +  [w°(Y°  -  H)]^  =  0. 

Equations  (2.10)  are  the  equations  of  the  nonlinear 
shallow  water  theory  for  three-dimensional  motion  over  an 
arbitrary  bottom.   The  functions  u  ,  w  ,  and  Y  depend 
upon  oc,    Y>  andr,  and  the  additional  condition  u  ^  =  w^ 

y        OS 


(from  equations  (2.4))  must  be  satisfied.   The  vertical 

veloci 

(2.9). 


velocity  v°  is  zero,  and  rr°  is  determined  by  equation 


3.   The  linear  shallow  water  theory 

The  linear  shallow  water  theory  csin  be  obtained  from 
equations  (2.10)  by  assuming  that  Y  =  a  +  £,  where  a  is 
constant  and  6,  u,  w  are  small.   If  quadratic  terms  in 
these  quantities  are  omitted,  equations  (2.10)  yield 


(3.1)      e^^  =  [e^(a  -  H)]^  +  [6^(a  -  H)]^- 

This  is  the  eqviation  for  the  surface  elevation  in  the 
linear  shallow  water  theory.   The  velocities  can  be  ob- 


tained  from  the  equations  u°  =  B.^    and  w  =  6   and  the 

o 

pressvire  from  equation  (2.9).   The  vertical  velocity  v 

is  zero. 

For  two-dimensional  motion  (e.g.  in  a  rectangular 
canal  of  constant  width)  w°  is  zero  and  therefore  u  and 
Y°  are  independent  of  y.   The  nonlinear  shallow  water 
theory,  equations  (2.10),  then  becomes 


12 


(3.2)  u°  +  u°u^  +  Y°  =  0, 

Y°  +  [u°(Y°  -  H)]^  =  0, 

C/  06 

with  V  still  zero  and  -rr^  given  by  equation  (2.9).  These 
are  the  well-known  differential  equations  of  the  nonlinear 
shallow  water  theory,  for  two-dimensional  motion,  first 
obtained  in  this  way  by  K.  0.  Priedrichs  [9].   The  linear 
shallow -water  theory  equation  (3.1)  becomes,  for  two- 
dimensional  motion. 


(3.3)  £^^  =  [e^(a  -  H)]^. 

4.   Higher  approximation^ 

In  order  to  get  more  accurate  solutions  than  those 
yielded  by  the  ordinary  shallow  water  theory,  we  continue 

the  development  in  powers  of  d .      This  is  done  in  the  pres- 

2      3 

ent  section  for  terms  in  of  and  d    ,   but  only  for  the  case 

of  two-dimensional  motion.  First,  considering  terms  in 
cf  in  eqxiations  (2.3)  and  making  use  of  the  results  ob- 
tained above  for  the  first  approximation,  one  has 

(4.1)  ^l  =  -v|, 

XT    +    U°U       +    U    U°     +    W       =    -V    U- , 

Z  00  oc  oc  p 

V^    +   U°V      +   TTq    =    -V    v„  , 

P         ot 

Y^  +  u°r^   +  uV   -  vV  =  v^  at  p   =  Y°, 

r  a  a,         p 


TT-"-  +  Tr°Y^  =  0  at  p    =  Y°, 


u  H     =  V  at  B   =  H. 
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The  third  and  fourth  equations  above  can  be  inte- 
grated if  equation  (2.8)  is  used  for  v  .   One  finds 

(4.2)  rrl=pV2(u°^.u%°^-uf) 

-^  P[^°(^°H)^  -  u°(u°H)^^  -  (u°H)^]  +  g(a,er), 

(4.3)  u^  =  -pVs  u°^  +  P(U°H)^^  +  f(cfc,^;). 

The  functions  f{oc,z)    and  g(a,?;)   are   so  far  \indetennined. 
Using  eqxaation    (4.3)    in   the  first  of   equations    (4.1) 
gives 

(A   A^          v2   -    (P^   -  h'^)    ,^o            ii!_-j£lf    OhI 
(4-4)  ^     -  -^ — 6 "^ococcc-  -^ — S ^"^  ^Vcfcct 

-    (p    -  H)f^    +  v}{<x.,E,z)E^. 

The  integration  constant  in  (4.3)  has  been  determined 
from  the  last  of  equations  (4.1) . 

Prom  the  sixth  of  equations  (4.1),  since  rr  =  -1, 
one  has 

(4.5)  TT-"-  =  Y-^   at  p  =  y°. 

This  equation  yields  g{oc,Z)    in  terms  of  i     as 

o2  2 

-1   Y   /  o     o  o    ,  o  \ 

(4.6)  g(a,r)  =  r"  -  V^'^o.t;  +  ^  ^aa  "  V  ^ 


-^°[->°n)^  --  (uH)^^-  (u^H)J. 

When  equations  (2.8),  (4.2),  and  (4.3)  are  used  in  the 
second  and  fifth  of  equations  (4.1),  the  following  equa- 
tions for  f (a,^)  and  Y^{oo,z)   sire   obtained: 
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(4.7)  f^  +  n\    +  u°f  =  -s^  -  ^[(u°H)f]^, 

1      O  1      O  1  _  Y°   -H^   O      Y°       -^     (Ojj)  /yO     . 

+  fV  ^c  -  Y°(^°H)^,  -  f](Y°  -H^). 

The  term  g^j^  is  to  be  eliminated  by  using  equation  (4.5). 
When  f  and  Y"^  are  found  from  equations  (4.7),  equation  (4.6) 
yields  g{oc,z)    and  equations  (4.2)  and  (4.3)  then  give  u  and  tt  . 
Equation  (2.8)  already  gives  v  .   Equations  (4.7)  are  the  equa- 
tions of  the  first  correction  to  the  shallow  water  theory  and 
will  be  called  the  equations  of  the  second  approximation,  the 
shallow  water  theory  being  the  first  approximation. 

Now  ret-urning  to  equations  (2.3)  and  equating  to  zero  tho 
coefficients  of  cf   ,    one  finds  for  the  third  approximation 

(4.8)  ul   =  -v^, 

2    2 

v^  .  VpV  .  VpV  =  u°Yf.u\^.u2Y^yf.UpVY°   at  p  =  Y°, 


-Y^  +  TT^Y-^  +  TT^  =  0  •  at  p  =  Y°, 


r2 


u^H^  =  v^  at  (3  =  H. 


Equations  (4.8)  for  the  third  approximation  couild  be  treated 
in  the  same  way  as  equations  (4.1),  but  this  will  not  be  done 
here . 
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5 .   Steady  state  solutions  and  waves  of  permanent  shape 

The  equations  for  the  first  three  approximations  in 
the  case  of  two-dimensional  motion  have  been  obtained. 
The  steady  state  solutions  of  these  equations  will  now  te 
investigated.   Such  solutions  yield  progressive  waves  of 
permanent  shape  when  they  are  referred  to  a  moving  coor- 
dinate system. 

Equations  (3.2)  of  the  first  approximation  (or  non- 
linear shallow  water  theory)  can  be  integrated  at  once  in 
the  steady  state.   When  integrated  they  are 

o 
(5.1)  ^  u°   +  Y°  =  e, 

u°(Y°  -  H)  =  m. 


where  m  and  e  are  constants. 

Inserting  the  expression  for  u  obtained  from  the 
second  equation  into  the  first  gives 

(5.2)  -2 +  Y°  =  e. 

2(Y°  -  H)'^ 

This  Is  a  cubic  equation  for  Y°  in  terms  of  H,  and 
shows  that  if  H  is  a  constant  (i.e.  the  bottom  is  hori- 
zontal) then  the  only  continuous  steady  state  is  Y  = 
constant  and  thus  u°  =  constant.   Since  the  equation  has 
one  root  for  which  Y°  -  H  is  negative,  and  since  this  is 
physically  impossible,  there  are  only  two  possible  con- 
tinuous steady  states  for  given  values  of  m  and  e.   There 
are  also  discontinuous  steady  state  solutions  consisting 
of  these  two  constant  solutions  pieced  together." 


""'  There  may  be  many  points  of  discontinuity,  and  at  each 
such  point  the  solution  changes  from  one  of  the  two 
constant  states  to  the  other.   However,  if  it  is 
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The  steady  states  given  by  the  first  approximation 
or  shallow  water  theory  having  been  obtained  for  the  case 
of  a  horizontal  bottom,  we  proceed  to  investigate  the 
second  approximation  to  these  solutions.   Equations  (4.7) 
become  in  this  case  (from  equation  (4.3),  f  =  u  ) 

(5.3)  ^°^i+Y^=0> 

(Y°  -  H)u^  +  u°y^  =0. 

The  only  solution  of  these  homogeneous  linear  equa- 
tions is  u^  =  Y^  =  0  and  therefore  u  =  constant, 
1     =  constant,  unless  the  determinant  of  the  coefficients 
is  zero.   The  determinant  is  zero  if  (u°)   =  Y°  -  H,  and 
in  this  case  Y~  =  -u  u^^  where  u^  is  arbitrary.   This  con- 
dition means  that  u  =  /g«deptW,  which  is  the  critical 
speed  or  propagation  speed  in  the  shallow  water  theory. 
Integrating  and  using  equations  (4.2)  and  (4.6)  for  tr 
and  equation  (2.8)  for  v  gives  the  solution 


(5.4)   v^  =  0,   tr^  =  r"-  =  -u°u^  +  c   if  (u°)^'=  Y° 


H 


where  u  is  so  far  an  arbitrary  f -unction  of  a  and  c  is  a 
constant. 

The  arbitrary  function  which  appears  in  the  second 
approximation  can  be  determined  by  considering  the  third 
approximation.   (This  circumstance  seems  to  be  rather 
typical  of  pertiirbation  solutions  of  boundary  value  prob- 
lems. )   Thus  the  steady  state  solutions  of  equations 


ass^umed  that  the  energy  of  a  small  element  of  fluid 
does  not  increase  on  crossing  the  discontinuity,  then 
only  one  discontinuous  solution  with  a  single  point  of 
discontinuity  is  possible.   For  this  solution,  the  flow 
is  from  the  shallower  to  the  deeper  side. 
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(4.8)  with  H,  u  ,  and  Y  constant  axid  the  second  approxi- 
mation given  by  equations  (5.4)  must  be  investigated. 
Prom  equation  (4.4)  one  has 

(5.5)  v^  =  -((3  -  H)u2^. 

Using  this  in  the  second  and  fouirth  of  equations  (4.8) 
gives 

(5.6)  TT^  =  -"°'^g-  "'  ^l^^^M. 

2       oca,        -^  \^/  ' 

Tl"-o  sixth  of  equations  (4.3)  then  yields 

(5.7)  Y^  =  Tr2(a,Y°)  =  Z2il(Y!^lJLL  u^^+  s(*). 

If  eqviations  (5.6)  are  used  in  the  second  of  equa- 
tions (4.8),  the  resulting  equation  may  be  integrated 
giving 

(5.8)  s  +  u°r  =  -  •|(u-'-)^  +  b 

where  b  5.s  con<?tant.   Similarly,  one  obtains  from  the 
third  of  equations  (4.8) 

(5.9)  v^  =  ]^«i»a-(p  -  H)^  -  r^((3  -  H) 


T/here  the  integration  constant  has  been  determined  from 
the  last  of  equations  (4.8).   liVhen  equation  (5.9)  is  used 
in  the  fifth  of  equations  (4.8),  the  resulting  equation 
may  be  integrated,  yielding 

(5.10)   uV  +  uV  =  -(Y°  -  H)r  +  (Y°  -  H)^  ^^^^^   +  j 
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where  j  is  a  constant. 

Equations  (5.7),  (5.8),  and  (5.10)  are  Inhomogeneous 
linear  algebraic  equations  for  the  determination  of  Y^, 
r,  and  s.   However,  the  determinant  of  the  coefficients 
vanishes  (because  (u°)   =  Y°  -  H)  and  therefore  the  terms 
independent  of  the  lanlcnowns  must  satisfy  a  condition  in 
order  that  the  equations  be  consistent.   The  condition 
may  be  obtained  by  multiplying  equations  (5.7)  and  (5.8) 
by  -u  and  then  adding  equations  (5.7),  (5.8),  and  (5.10). 
The  result,  after  Y"^  is  eliminated  by  use  of  equation 
(5.4),  is 

,_  T_.   -(Y°  -  H)"^   1   .  3   o,  lv2     1    ,.    o,v    ^ 

(5.11)  -^ g  '      u^^  +  ^  u  (u  )   -  cu  +  (j  -  u  b)  =  0. 

This  condition  determines  u  (a) ,  which  was  not  determined 

by  the  equations  of  the  second  approximation.  Y/hen  u  (a) 

is  found  the  solution  of  the  second  approximation  will  be 

complete,  and  equations  (5.7),  (5.8),  and  (5.9)  for  the 

third  approximation  can  be  solved.  However,  an  arbitrary 

fixnction  will  appear  in  this  solution,  and  it  would  be 

necessary  to  go  to  the  fourth  approximation  to  determine 

it. 

Equation  (5.11)  is  of  the  same  type  as  that  obtained 

by  Rayleigh  and  Korteweg  and  DeVries,  and  may  be  inte- 

6ul 

grated  if  it  is  multiplied  by  — — »  ,  and  the  result  is 

(Y°  -  H)^ 

ic,    ^o^   /„1^2  _    3u°   /„1n3      3c  ,\.2   ^   6(  j-u°b)  ,  1  .  „ 

(5.12)  (u  )   =  — ,r(u  ) t(u  )   +  — i-^ k   u  +q 

"^      (Y°-H)^       (Y°-H)^        (Y°-H)^ 

where  q  is  a  constant.   This  equation  may  be  integrated  in 
terms  of  elliptic  functions.   To  this  end,  let  the  roots  of 
the  cubic  on  the  right  be  £,   h, ,  and  hg  where  /  >  h^  >  hg 
if  all  roots  are  real.   If  the  roots  are  not  all  real, 
equation  (5.12)  has  no  non-constant  solutions  which  are 
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bounded  in  the  region  -co  <  oc  <    +00 .   Since  the  svira  of 

1  p 
the  roots  is  the  negative  of  the  coefficient  of  (u  )   di- 

1  3 
vided  by  that  of  (u  )  ,  the  roots  must  satisfy  the  condi- 
tion ^ 


(5.15)  /+  h^  +  hg  =  -^  . 

u 

Now  introduce  the  new  variable  X   ^7 

(5.14)  u^  =  h^  +  (hg  -  h^)    cos^  X 
Then  equation  (5.12)  becomes 

(5.15)  ^  =  A-^^1  -  k^  sin^  X 
where 


Integrating  equation  (5.15), 

(5.17)  a  -  a^  =  A  r   (1  -  k^  sln^  X)~'^^^   dX  • 

o 

The  origin  may  be  chosen  so  tlaat  oc     =  %  =  0,  which 
implies  that  the  origin  occurs  where  u-^  =  hp.   Thus 
a  =  LFiXt'k.)    where  P(X,k)  is  the  elliptic  integral  of  the 
first  kind  of  modulus  k.   The  inverse  function  is 
X  -   arc  cos  en  —(mod  k)  and  thus 

(5.18)  u-^  =  h^  -  (h^  -  h2)cn^  |(mod  k) . 

This  solution  is  periodic  in  a  if  /,  h.^,    and  hg  are 
real,  and  the  wave  length  is  then 
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(5.19)     X  =  2t     f  (1  -  k^  sln^  X)~'^^^   dX 


rTT 


=  2LF{^  ,  k)  =  2AP-L(k). 

If  the  roots  are  not  all  real,  the  solution  becomes 
infinite  as  a  approaches  infinity.  For  this  reason  only 
the  case  of  real  roots  will  be  discussed. 

From  equations  (5.4)  and  (5,18)  one  has 

(5.20)   TT"^  =  Y^  =  c  -  n^h^   +   u°(h^  -  h2)cn^  |(mod  k) . 

The  constant  c  may  be  determined  if  it  is  required 
that  the  mass  of  fluid  in  a  wave  length  equal  ?l(Y  -  H), 
the  amount  of  fluid  per  distance  %   ir  -^hc  first  app'"'""^'' - 

mat  ion.   Then  /   Y"^d(t  =  0,  which  leads  to 


/' 


'^  E  (k) 

(5.21)  c  =  u7  -  u°(/  -  hg)!^^  =  u°(/  +  h3_  +  hg). 

Here  E-,  (k)  is  the  complete  elliptic  integral  of  the  sec- 
ond kind,  and  the  last  equality  is  a  consequence  of  equa- 
tion (5.13) . 

The  second  approximation  has  now  been  determined  and 
it  has  led  to  a  two-parameter  family  of  solutions.  Now, 
combining  the  first  ancL  secoi*^  ttPpro::linatio*xi> ,  ono  laas 

(5.22)  V  =  v°  +  cfv^  =   0, 

u  =  u°  +  cTu"  =  u°  +  dh-j_  -  (<?h^  -  c(h2)cn  ^(mod  k) , 

Y  =  Y°   +  CJY^  =  Y°   +  u°(cfi  +  cfhg) 

+  u°(c?h^  -   cJh2)cn^  ^(mod  k) , 

TT  =  Tr°   +  <?Tr^  =  Y°   -  p,  +  u°(di  +  dhg) 

+  u°(cJh,   -  rfh„)cn^^(mod  k) . 

1  <5  fl 
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2  oc 
Since   en     —  varies   betv/een  zero  and  one,    the  maximum 

and  minimum  values   of  Y  are   given  by 
(5.23)  Y^^^  =  Y°   +  u°(cj/  +  ah^), 

Ymin=Y°   +u°(cr/-fc.h2). 

These  equations  v/ill  be  used  to  eliminate  cJh,  and 
rfhp  from  equations  (5.22).   Also  it  will  be  assumed  that 
H  =  0,  which  means  that  vertical  distances  are  measured 
from  the  bottom,  and  that  Y  =1,  which  Implies  that  h  is 
the  average  height  of  the  surface.   Then,  upon  reintro- 
ducing the  original  coordinates,  equations  (5.22),  (5.21), 
and  (5.19)  become  (the  substitution  L  =  cf/  has  been  used) 


(5.24)    V  =  0, 

u  _  '^max   -.    .'^max   '?minv„„2rX  ^-ci  nr\T  r^^A   v^ 
-  — ^ L  -  (      ^ E~'    '•X  2P-j_(k)]  (mod  k) , 


n  =  ^min  +  (^max  "  ^min^^^^tf  2P^(k)]  , 
|=Ap^(k)(2L.l-^)"'^', 


where 


^mnx   ^min 
2L  +  1  >  22g5  >  ?2l£  ,   0  <  k^  = '^  ""?"  <  1, 


and 


(2L  +  1  -  ^)E3^(k)  =  (2L  -  2  -  l?2g5  .  !^)F3_(k). 

Equations  (5.24)  represent  the  steady  state  solutions 
as  given  by  the  first  and  second  approximations.   These 
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solutions  have  already  been  discussed  In  the  Introduction. 
V/hen  A  =  CO ,  equations  (5.24)  yield  the  solitary  waves 
given  by 


(5.25)    v  =  0,     X  =  oo, 


T=i  -  T5( — Tr-+l)-( — TT--!)  sech  -^{ — r— -1) 


^  -  1  +  (-^  -  1)  sech  -2H(-Tr  -  ^)    ' 


1-  =  -n  -  J. 
PS    '   ^ 

This  solution  has  also  been  discussed  in  the  introduction. 
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